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We present a theory of the linear and nonlinear optical characteristics of the insulating phase 
of the Falicov-Kimball model within the self-consistent mean-field approximation. The Coulomb 
attraction between the itinerant d-electrons and the localized /-holes gives rise to a built-in coherence 
between the d- and /-states, which breaks the inversion symmetry of the underlying crystal, leading 
to: (1) electronic ferroelectricity, (2) ferroelectric resonance, and (3) a nonvanishing susceptibility 
for second-harmonic generation. As experimental tests of such a built-in coherence in mixed-valent 
compounds we propose measurements of the static dielectric constant, the microwave absorption 
• spectrum, and the dynamic second-order susceptibility. 
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C/3 ■ I. INTRODUCTION 

In this paper, we present a theory of ferroelectricity originating from an electronic phase transition, in contrast to 
the conventional displacive ferroelectricity due to a lattice distortion |IJ. The electronic ferroelectricity occurs in a 
strongly correlated electron system, namely, the insulating phase of the Falicov-Kimball model. 

The Falicov-Kimball (FK) model was introduced originally Q as a simple model to explain the metal-insulator 
transitions observed in certain transition-metal and rare-earth compounds. The model has since been applied exten- 
sively to the mixed-valent compounds and heavy-fermion materials. The FK model introduces two types of electrons: 
■ itinerant d-electrons and localized /-electrons. The valence transition is driven by the on-site Coulomb repulsion 
between the d- and /-electrons. A d-f hybridization term may or may not be added to the model. The theoretical 
solutions for the ground state of the FK model can be divided into two classes. On the one hand, solutions with 
/-occupation as a good quantum number |3|,|j do not have a built-in coherence between ci-electrons and /-holes. On 
the other hand, solutions such as the self-consistent mean-field (SCMF) solution j|] and the electronic polaron || do 
have a built-in coherence between d-electrons and /-holes. 

The built-in coherence of the SCMF solution breaks the inversion symmetry of the FK Hamiltonian in the following 
way. The many-electron eigenstates of the Hamiltonian can be classified into even and odd parity states. The SCMF 
ground states with even and odd parities are degenerate. A linear combination of the even and odd parity states 
forms the appropriate ground state for the limit of a vanishing electric field. We shall argue that such degenerate 
O ■ ground states can exist for solutions with /-occupation as a good quantum number. Thus, the inversion symmetry 
breaking is not limited to the SCMF solution. 

The primary purpose of this paper Q] is to give a detailed account of the linear and nonlinear optical characteristics 
of the SCMF solution. The inversion-symmetry-broken ground state possesses the following distinctive properties: 
(1) electronic ferroelectricity, (2) ferroelectric resonance, and (3) a nonvanishing susceptibility for second-harmonic 
generation. Solutions of the model without built-in coherence do not have these properties. As experimental tests 
to distinguish between the two classes of solutions we propose measurements of the static dielectric constant, the 
microwave absorption spectrum, and the dynamic second-order susceptibility of a mixed-valent compound, for example 
SmB6. 

In recent years, Four- Wave-Mixing (FWM) spectroscopy has emerged as a powerful tool for studying coherence 
in optically pumped semiconductor systems |&-Cd]. In a three-beam FWM experiment, two incoming beams of 
wavevectors ki and k2 set up a transient polarization grating. A third incoming beam of wavevector ka diffracts off 
this grating to produce an outgoing signal in the direction ki = k3 + k2 — ki. Being a third-order process, FWM 
is allowed in media with or without inversion symmetry. We pose the question: what happens if the system being 
probed already has a polarization built into it by nature? An example of such a system is the SCMF solution of the 
FK model resulting in the Bose-Einstein condensation of d-f excitons. 
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As shown below, the built-in polarization leads to a nonlinear optical response to second order in the external field. 
The mixed-valent system has a nonvanishing susceptibility (2u>, u>, uS) for second-harmonic generation. The built- 
in polarization replaces one of the incoming beams of the three-beam FWM experiment. In crystals with inversion 
symmetry, second-harmonic generation is forbidden in the electric-dipole approximation. In the mixed-valent system 
the built-in polarization breaks the inversion symmetry, allowing second-harmonic generation to take place. We present 
a calculation of the second-harmonic susceptibility of a model mixed-valent system within the SCMF approximation. 
The second-harmonic susceptibility is directly proportional to the built-in coherence A, showing that second-harmonic 
generation can be used as a test of d-f exciton condensation in mixed-valent compounds. 

The existence of a built-in polarization in the ground state also means that, according to the SCMF theory, mixed- 
valent compounds are ferroelectric. Whereas in commonly known ferroelectrics the built-in polarization is due to 
the relative displacement of positive and negative ions, the ferroelectricity in mixed-valent compounds is of purely 
electronic origin. Apart from possible Jahn- Teller distortions as a result of the electronic polarization, the valence 
transition does not involve a change in the crystal structure. As shown below, the valence transition is accompanied 
by a divergence of the static dielectric constant at the critical value of the /-level energy. The divergence of the static 
dielectric constant should be observable in real mixed-valent compounds, for example by varying the external pressure 
or temperature. 

The continuous symmetry associated with the phase of A leads to a Goldstone mode in the excitation spectrum of the 
mixed-valent compound. In the pseudo-spin picture, the Goldstone mode corresponds to a uniform precession of the 
pseudo-spins around the z-axis. Ferroelectric resonance occurs when an ac electric field is applied whose frequency 
coincides with that of the Goldstone mode. This phenomenon is the electric analogue of magnetic resonance in 
ferromagnetic insulators. The ferroelectric resonance frequency is proportional to the square root of the effective bias 
field, and may depend on the sample shape, domain structure, and crystal fields. In real mixed-valent compounds, 
ferroelectric resonance should occur in the microwave regime. 

The remainder of the paper is organized as follows. In section |H| we discuss the SCMF solution for the ground state of 
the FK model. We show that the SCMF ground state has a built-in polarization which breaks the inversion symmetry 
of the FK Hamiltonian. We calculate the static dielectric constant of a model mixed-valent system using mean-field 



theory. In section III we calculate the linear susceptibility of the model system using the pseudo-spin formalism. 



We obtain an analytical expression for the linear susceptibility in a uniform static electric field. We determine the 



ferroelectric resonance frequency and analyze the shape of the infrared absorption spectrum. In section IV we compute 
the second-harmonic susceptibility of the model system. In section [v| we compare the results of the model calculation 



to experimental data and propose experimental tests of coherence in real mixed-valent compounds. In section VI we 



discuss the possibility of ferroelectricity, ferroelectric resonance, and second-harmonic generation for solutions of the 



FK model other than the SCMF solution. The main results are summarized in section VII 



II. MODEL 



Ignoring the electron spin, the FK Hamiltonian is 

H = E £kd lA + E' f fth + E + h-c. 



k k k 

k,k',q 

Here d^ creates a rf-electron of momentum k and energy e^, and fl creates an /-electron of momentum k and energy 
E'f. The parameter U is the Coulomb repulsion between the d- and /-electrons, Vk is the hybridization energy, and N 
is the number of sites. We assume that the d-band and the /-level are derived from d- and /-orbitals on the same site. 
We have chosen the d- and /-orbitals to yield a finite dipole moment between them in the ^-direction. For simplicity, 
we consider a model system with a d-band of bandwidth W and constant density of states po = 1/(2W). 



A. SCMF solution 



The SCMF solution is analogous to the BCS theory of superconductivity except that the pairing now occurs between 
a d-electron of momentum k and an /-hole of momentum — k (Ref. [pT[ ). Whereas a Cooper pair carries charge, an 
electron-hole pair is neutral. The SCMF solution therefore describes an insulator, rather than a superconductor. 
Pairing between electrons and holes may also occur in a semiconductor placed in an intense coherent laser field jL2) . 
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In that case, the pairing is externally enforced by the pump field. The mean-field decoupling of the Hamiltonian ([[]) 
yields the effective one-particle Hamiltonian 



H = ]T(e k + Un f )d f k d k + (E' f + Un d ) ^ /+/ k 

k k 

+ £(Vk - A)4/k + h.c. - NUn f n d + ^|A| 2 , (2) 



U' 

K 

where A = ^ £k(^l/k d kM is the gap parameter, and n d = ± £ k (V'l4AlV') and n/ = ^ Ek(^l/k/k|^> = 1 
are the e?-band and /-level occupancies respectively. The SCMF ground state is \tp) = Yik( u k + u kdt/k)|0), where 
|0) is the state with no /-holes (the normal state), and Uk = cos («k — sin h&k) is the probability amplitude for 
the pair state (k, — k) to be occupied (unoccupied). The gap parameter A and the /-level occupancy rif must be 
determined self-consistently from 

A=££— , (3) 
2E k ' 

k 



n f 



k v k 



Here E^ = + |A — Vk| 2 is the quasiparticle excitation energy, with £k = ^(^k ~ Ef), where Ef = E'f + Un is the 
renormalized /-level energy, and n = n d — nf is the inversion. Since neither Ef nor are known from first principles, 
we shall treat Ef as the materials parameter. Eqs. @ and (^) are Eqs. (11) and (10) of Ref. || at temperature T = 0, 
with a k-dependent hybridization. If the crystal has inversion symmetry, the hybridization must satisfy = — Vk- 
The assumption of a k-independent hybridization in, among others, Ref. is therefore incorrect. If, instead, we 
assume nearest-neighbour hybridization, we find that Vk is odd in k and purely imaginary. As can be seen from 
Eq. (||) , the imaginary part of A then vanishes due to the cancellation of terms with ±k. The real part of A is given 
by the BCS gap equation 

N ^ 2E U ' y ' 

k 

where E^ = + | | 2 + A 2 . Calculation shows that a sufficiently strong hybridization can destroy the gap. In the 
following we consider the limit where Vk is negligible compared to U. 
The solution of Eqs. (§) and (@) 

for our model system is shown in Fig. |l|. The figure shows the gap parameter A 
and the /-level occupancy n/ as a function of the /-level energy Ef. The gap parameter A is the order parameter of 
the valence transition. When the /-level is far below the bottom of the d-band, the system is in the normal state with 
no /-holes, and A = 0. As the /-level is moved upward past the critical value Ef = —E c , where E c = W coth(W/U) 
(in a real material this is achieved by applying pressure or by alloying) , A becomes nonzero and the system undergoes 
a valence transition. In the mixed- valent state, the /-level occupancy rif lies between and 1. The gap parameter 
reaches the maximum value A(0) = W/[2 sinh(W/U)] when the /-level lies at the center of the <i-band (Ef = 0). 
This is the half-filling point nf(0) = 1/2. For still higher Ef, the /-level gradually empties out into the <i-band. The 
solutions transform according to A(Ef) = A(—Ef) and rif(Ef) = 1 — n(—Et). At Ef = E c , the system returns to 
a normal state with no /-electrons. Since nf(Ef) has no discontinuities, the valence transition is continuous for all 
values of the Coulomb repulsion U. 



B. Inversion symmetry breaking 

The key feature of the SCMF ground state is that it breaks the inversion symmetry of the FK Hamiltonian. 
In general, symmetry breaking occurs when the ground state has a lower symmetry than the Hamiltonian. As a 
well-known example in another area, the Ising model -fusing = ~JJ2<ij> SfSj is invariant under Sf — > —S*, all i. 
Yet the ground state is either one of the broken-symmetry states with built-in magnetization ±M Z . The sign of the 
magnetization is selected by applying an infinitesimal bias field H z , and then choosing the lower-energy state. 

In the case of the SCMF solution of the FK model, the inversion symmetry in a d-f site is spontaneously broken 
due to the pairing of ci-states of even parity with /-states of odd parity. Applying the inversion J to \ip) gives the 
state 
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j\^ = H(-K+v* k dthm, 

k 



(6) 



which is linearly independent of l^). The form of the inversion image follows from the even parity of the (i-orbital 
and the odd parity of the /-orbital and M_k = u£ and v_k = The states \ip) and J\if)} have built-in polarizations 
(il)\P\ip) = P (0) and {ip\JPJ\ip) = -P(°), where 

k 

is the polarization operator, fi is the interband dipole matrix element, which for simplicity we take to be independent 
of k, and ft is the volume. One can show that in the thermodynamic limit N — > oo, and J\tp) are orthogonal and 
uncoupled by H . The proof relies on the fact that the infinite product Ilkd^kl 2- l u k| 2 ) is zero, since |(|t>k| 2 — l w k| 2 )| < 1 
for almost all k. 

Because H is invariant under inversion, and J\ip) are degenerate. As in the case of the Ising model, the correct 
ground state is selected by lifting the degeneracy with a bias field E, and then choosing the lower-energy state. This 
yields the state \ip) with built-in polarization in the direction of E. We call this direction the z-direction. (Without 
crystal-field terms, the z-direction has no definite orientation with respect to the crystal axes.) Since fj, z is real, 

Pz QS> = N(i z (A + A*)/(ft£7), where A is the built-in coherence. The built-in polarization vanishes in the normal state 
where A = 0. 



C. Electronic Ferroelectricity 

In commonly known ferroelectrics such as BaTiC>3, the ferroelectric transition involves a change in the crystal 
structure. In the ferroelectric phase, the positive ions are displaced relative to the negative ions, leading to a permanent 
electric dipole moment. The displacive ferroelectric transition occurs when the transverse optical (TO) phonon 
frequency vanishes at some point in the Brillouin zone. In an electronic ferroelectric, the ferroelectric transition 
involves a change in the electronic structure rather than the crystal structure. (Here we neglect the electron-phonon 
coupling, which may cause a lattice distortion as a secondary effect of the transition.) Instead of a vanishing of the 
TO phonon frequency, the d-f exciton energy goes to zero at the critical value of the /-level energy. The built-in 
polarization of an electronic ferroelectric is of the order of 10 /LtC/cm 2 , comparable to the built-in polarization of 
perovskites Q). 

Since the built-in polarization is continuous at Ef = ±E C (see inset of Fig. ^|), the valence transition is a second- 
order ferroelectric transition. In general, second-order ferroelectric transitions are accompanied by a divergence of the 
static dielectric constant in the direction of the spontaneous polarization. For a temperature-driven transition, the 
dielectric constant diverges as (T — T v )~~> above T r , and as (T r — T)~"' below. Almost all known ferroelectrics have 
7 = 1 (Curie- Weiss law). Observed values of 7' range from 1/2 in TGS to 1/6 in SbSI |T^| . 

Here we calculate the static dielectric constant of the model mixed- valent system at temperature T = 0, us- 
ing mean-field theory. The dielectric constant in the z-direction is given by e Z2 = 1 + ^x'zz , where Xzz = 
lhriB^o dPz/dE z is the static susceptibility. For A real, the polarization is P z ^ = 2N/j, z A/(QU), so that 
Xzz = [2 A/i z /(£![/)] lim£ z _>o dA/dE z . The bias field E z leads to an additional term ~^ z E z ^2^d\ L f k + h.c. in 
the effective one-particle Hamiltonian. The self-consistency equation for A then becomes 

A = ^V A + ^, (8) 
N 2E k y ' 

k 

where Ek = + (A + [i z E z ) 2 . The susceptibility is obtained by implicit differentiation of Eqs. (§) and (@) with 
respect to E z . The susceptibility of the normal state (A = 0) is 

v (i) = 2N(il arccothflfl/l/WQ 

Xzz ft W- Uarccoth(\E f \/Wy { ) 

and the susceptibility of the mixed- valent state (A > 0) is 

x « = --M U + l + ^A(0) \ 

Xzz flU \ 4A 2 A(0)[1 +4A 2 A(0)] + 4A 2 B 2 (0)y 1 ' 
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Here A(0) and -B(O) are given by Eqs. ( |5^ ) and ( |53| ) of the Appendix. Fig. || shows the dielectric constant of the 
model system as a function of the /-level energy. The dielectric constant diverges at \EA — E c . From Eq. @ we 

find x*z oc (\Ef \ — -Ec) -1 as \Ef \ approaches E c from above, and from Eq. ( [To| ) we find xi* (-^c — l-^/l) -1 as l-^/l 
approaches E c from below. Thus, the critical exponents according to mean- field theory are 7 = 7' = !. 



III. LINEAR OPTICAL RESPONSE 



We first consider the linear optical response of the mixed- valent system. The SCMF solution predicts an energy gap 
2A in the absorption spectrum. The gap is 2A because the incoming photon must create two quasiparticles, just as in 
a superconductor. Far-infrared transmission and reflectivity spectra fll|-|l6f|, as well as electron tunneling spectra [fl7|| , 
show energy gaps of several meV in a number of mixed-valent compounds. The crucial difference between the mixed- 
valent compound and the superconductor is this: in the superconductor, the pairing occurs between two electrons, 
whereas in the mixed-valent compound the pairing occurs between an electron and a hole. This has important 
consequences for the coherence factors that enter the response of both systems to different external probes. For 
example, it is a textbook result ]18[ ] that a clean superconductor at temperature T = cannot absorb electromagnetic 
radiation because the coherence factor M p u p _|_ q — w p u p + q vanishes for zero photon momentum q. For the mixed-valent 
compound the coherence factor entering the electromagnetic absorption is w p it p + q — u p w p + q , which remains finite for 
zero photon momentum. The coherence factor entering the electromagnetic absorption of the mixed-valent compound 
is the same as the coherence factor entering the acoustic attenuation rate of the superconductor JToj ] . 

We calculate the linear response of the mixed-valent system to an ac electromagnetic field in the presence of a dc 
bias field. The bias field serves to select the direction of the built-in polarization. In a real material, the bias field 
is provided by the crystal field or the depolarization field due to the sample boundary. We treat the interaction of 
the mixed-valent system with the ac electromagnetic field in the electric-dipole approximation. The interaction term 
in the Hamiltonian is ifi n t = -fc^Ek4A + h.c, where £ z is the component of the ac electric field along the 
z-direction (i.e. the direction of the built-in polarization). Only the z-component of the ac electric field couples to 
the channel in which the pairing takes place. The optical signatures of d-f exciton condensation occur only in this 
channel. We ignore the response of the remaining optical channels. 

The ac electric field sets up a polarization P z in the material, which in general can be a complicated nonlinear 
function of E z . In linear response, we expand P z in powers of £ z and keep only the first-order term: pj 1 ^ = xiVf z . 
In the electric-dipole approximation, the linear susceptibility xiV depends on the photon frequency ui but not on 
the photon momentum q. The quantity measured in experiments is the reflectivity spectrum or the transmission 
spectrum. From these one can extract the optical conductivity o~ zz {iS) by Kramers-Kronig analysis. The optical 
conductivity is related to the linear susceptibility by a zz (uj) = —iu>Xzz(u)- 



A. Optical Bloch equations 



We have calculated the linear susceptibility \ zz both from the Kubo formula and from the optical Bloch equations. 
The pseudo-spin picture gives a nice physical description of the linear and nonlinear responses of the system as 
precessional modes of the pseudo-spin vector Sk = (S x ^S y ^S z ^). The optical Bloch equations describe the time 
evolution of the pseudo-spin vector under the action of the ac electric field £ z . The components of the pseudo-spin 
vector are the expectation values of the pseudo-spin operators 

0-a.k = 4/k + /k d k, (11) 
0-y M = -i(4/k - /K). ( 12 ) 

o zM = 4d k - fth (13) 

in the ground state \ip). The equations of motion for the components of Sk follow from the Heisenberg equations of 
motion for the pseudo-spin operators (h = 1, i = x, y, z), 

Oi,k = ~i [ci,k, H + H int ] . (14) 

Working out the commutators, we find that the right-hand side contains products of pseudo-spin operators <7i, k c i j,k / 
(i ^3). The products occur because of the Coulomb interaction term in the Hamiltonian (fil). A closed set of equations 
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is obtained by replacing the average of products (ip\ui^aj^'\ij^) by the product of averages {ip\<Ji^\'4 1 ) {4'\ <J jM'\i } ) — 
SikSj t k'. This gives the optical Bloch equations 

S k = (H k - M k ) x S k , (15) 

where H k = (— 2fi z (E z + £ z ), 0, e k — E'*) and M k = ^X)k^k- The symbol x represents the vector cross product. The 
optical Bloch equations describe the coupled motion of a collection of TV pseudo-spins. Each pseudo-spin precesses 
around a local "magnetic" field H k — Mk, which is the sum of an external field Hk and an average internal field 
— M k , where Mk is the pseudo-magnetization. 



B. Stationary solution 

In the absence of the ac electric field the optical Bloch equations have a stationary solution S k . The stationary 



solution is obtained by setting S k °^ = in Eq. (15) 



0=(H k o) -M k o) )xS k o) . (16) 

Here H^ ' = {—2fj, z E z , 0, e k — E'*), and M k °^ = -j^ ^k°' ■ ^ n * ne stationary state each pseudo-spin is lined up with 
the local "magnetic" field. Then there are two possibilities: S k is either parallel or antiparallel to H^ - 1 — M k . The 

state with S k antiparallel to H k °^ — M k °^ has the lower energy. Thus in the ground state all pseudo-spins point in 
the direction opposite the local "magnetic" field. 

For zero bias field (E z = 0), Eq. ( pif ) is invariant under rotation about the z-axis. If S^ ' is a solution, then 
so is R z (<p)S^\ where R z {4>) is a rotation about the z-axis over an angle <fi. The angle <ft is the phase of the gap 
parameter A. For nonzero bias field, S k °^ lies in the x-z plane and the gap parameter is real. Introducing spherical 
polar coordinates, S k = (sin 6> k , 0, cos # k ) and M k °^ = (2A, 0, Un). The tilting angle k is the angle between S k °^ 

and the positive z-axis. The magnetization M k °' must be determined self-consistently from M^ " 1 = S ^ k S^ "*. The 
z-component of this equation gives Eq. (||), and the x-component gives Eq. (||). The stationary solution of the optical 
Bloch equations is the SCMF solution given in section [n| 



C. Pseudo-spin precession 

To calculate the linear susceptibility from the optical Bloch equations we expand the pseudo-spin vector and 
pseudo-magnetization to first order in E z : S k = S k 0) + and M k = M k 0) + M k 1} . Also, H k = H k 0) + H k 1} , with 

H k ^ = (— 2\x z £ Zl 0, 0). We substitute the expansions into the Bloch equations and collect terms of the same order in 
the ac electric field. To zeroth order we recover Eq. (tut). To first order we get 



S k 1} - (H k 0) - M k 0) ) x S k 1} - S k 0) x Mf) = H k x) x S k 0) . (17) 

Taking the inner product with the stationary solution gives S^ 1 "* • S k °^ = 0. The ac electric field causes the pseudo- 
spin vector to precess about the stationary direction. With S k tilted away from the z-axis, the precession involves 
variations in all three cartesian components of S k . The problem is simplified by working in the spherical polar 
coordinate system. In spherical polar coordinates, the stationary solution is the unit vector in the radial direction 
e r . The precession about e r is decomposed into components along the polar and azimuthal units vectors eg and e^: 
sj^ = Sg^eg + Si , k e,£. One must remember that e r and eg vary with the tilt angle # k , while e^ is fixed since S k °^ 
lies in the x-z plane. The equations of motion for Sg 1 ^ and are 

S$-21kS$ + M$ = (18) 
S$l + 2 £kSg - Afg = 2fi z £ z cos k . (19) 

Here = S^ k , and M^ k = jj ^ k , cos(0 k — &k')Sg^,. The appearance of the cosine factor in Mg 1 ^ is due 

to the variation of eg with # k . In deriving Eqs. ( |l8| ) and (19) from Eq. ( p~7| ) we have used the fact that 5^ k = 0, which 



follows from 

S*W = o an( j the initial condition S^l(t = 0) = 



G 



D. Noninteracting quasiparticles 



It is instructive to first calculate the absorption spectrum neglecting the Coulomb interaction between the optically 
excited quasiparticles. In the pseudo-spin picture, this corresponds to setting — Mg^ = in Eqs. ( |l8| ) and 

(|l9|). The linear susceptibility is \zz — P^ /£z, where P^ — (N^l z /Q,) ^ k Sg 1 ^ cos 6*k is the first-order polarization. 
Solving Eqs. @ and (|l|) for S$., we obtain 

(i) = _Mr ( U ^ k »k : k 
Xzz CI ^ \ uj - 2E k uj + 2£ k 

where rik,k' = Uk^k' — vuvw is a coherence factor. The frequency uj is understood to have a small positive imaginary 
part S. The physical origin of the coherence factor n k y can be understood as follows p0| : In the absorption process, 
an incoming photon of momentum k — k' creates a quasi-electron of momentum k and a quasi-hole of momentum 
— k'. This can be done in two different ways: (1) If the pair states (k, — k) and (k', — k') are initially empty, by adding 
a ci-electron in k and adding an /-hole in — k'. This process has amplitude u k u k /. (2) If the pair states (k, — k) and 
(k', — k') are initially occupied, by removing an /-hole from — k and removing a d-electron from k'. This process has 
amplitude — v^v^i. The overall amplitude for the creation process is Ukitk' — ^k^k' = ^k.k'- 

For the model system the absorption spectrum can be found analytically. When the /-level lies inside the (i-band 
(\Ef\ < W), the energy gap is 2A, and the absorption rises as \J uj — 2 A above threshold. When the /-level lies 
outside the d-band (\Ef \ > W), the gap is \/{\Ef\ — W) 2 + 4A 2 , and there is a discontinuous jump in the absorption 
at threshold. The single-quasiparticle result for the absorption spectrum at half-filling (Ef = 0) is shown by the 
dash-dotted line in Fig. [| 




E. Final-state interaction 



The quasiparticles created in the optical transition interact via the Coulomb potential U. In the normal state 
(A = 0) the final-state interaction leads to the Wannier exciton. In the mixed-valent state (A > 0) the final-state 
interaction leads to ferroelectric resonance and a threshold singularity in the infrared absorption spectrum. In the 
pseudo-spin picture, the final-state interaction corresponds to -Ml 1 ^ and Mg 1 ^.. For a separable Coulomb potential, 
Eqs. © and @ can be solved analytically. The pseudo-spin components are given by 



(21) 

rn , , n L^ COS 2 (^)r(k>) 

F(k, U ) = - cos 9 k - - £ — _ 

k' 

N ^—f u + 2E k , { 6) 

is a vertex function. In diagrammatic terms, the vertex function is the sum of all ladder diagrams contributing to the 
propagation of the quasiparticle pair. 

In the normal state (#k = 7r) Eq. (|2^) yields the Wannier exciton. The pole of the vertex function gives the exciton 
binding energy E^ = W coth(W/U) — W . The absorption spectrum of the normal state consists of an exciton line at 
uj = Eg — Eb, and a continuum between ui = E g and uj = E g + 2W, where E g — —W — Ef is the d-f band gap. The 
final-state interaction enhances the absorption at uj = E g . 

In the mixed-valent state the vertex function has three components: 

r(k,w) = -cos6» k r 1 (w) + sin6»kr 2 (w) + r3(w). (24) 



r(k, w ) r(k,- w )\ 



where 



The components Ti(u>) and T2(ui) are even in uj, while Ts(uj) is odd. Substitution of Eq. (24) into Eq. (23) yields a 
set of of three algebraic equations for the vertex components: 
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{uj 2 - M 2 )A{uj) + R MB{uj) ujB(uj) \ ( T x {u) \ (l 

MB(u) 1 + M 2 A(uj) ujMA(uj) T 2 {uj) = ) . (25) 

ujB(uj) ujMA{uj) uj 2 A(uj) + R J \T 3 (J) J \0 

Here M = 2 A + /i z E z and R — [i z E z j (A + fi z E z ). The functions A(uj) and B(ui) are given by 

A H = n E 2 £ k (c - 2e\)(uj + 2£ k ) ' (26) 

S H = jvE 2 £ k (co - 2E k )(u + 2E k ) ' (27) 

For the model system and can found analytically. The explicit expressions are given in the Appendix. 

Solving Eq. j25| ) for Ti(uj) gives 

r . , = (uj 2 + RM 2 )A(uj)+R 

U ; [(uj 2 -M 2 )A(uj) + R][(uj 2 + RM 2 )A(uj)+R}~(uj 2 +RM 2 )B 2 (uj)' 1 ' 



The linear susceptibility is 



^H = ^(riM-i). (29) 



F. Infrared absorption spectrum 

We first consider the absorption spectrum in the absence of a static electric field. We calculate the continuum 
absorption in zero bias field, setting M = 2 A and R = in Eq. (p8|). For an energy gap of several meV, the 
continuum lies in the far infrared. The linear susceptibility in zero bias field is 

'gw-^( ^-^L)-^) - 1 )- (30) 

The solid lines in Fig. || show the imaginary part of Eq. (^0|) for several values of Ef. When the /-level lies far 
below the bottom of the (i-band (top curve), the continuum absorption of the mixed- valent state is very similar to 
the continuum absorption of the normal state. There is a slight enhancement near the energy gap. As the /-level 
approaches the bottom of the d-band, the enhancement becomes more and more pronounced. When the /-level lies 
within the <i-band, the spectrum has a threshold singularity at uj = 2A. The spectrum for Ef above the center of 
the d-band is the same as the spectrum for — Ef below it. From Eq. j30| ) we find that when \Ef\ < W, the threshold 
singularity is e^ 1 ^ 2 9(e), and when \Ef \ = W the singularity is e -1 / 2 ln~ (e)6(e), where e = uj — 2A. When the /-level 
lies outside the <i-band the singularity is cut off because the energy gap is larger than 2A. The singularity is due 
to the final-state interaction between the optically excited quasiparticles. In the single-quasiparticle treatment, the 
absorption rises continuously from zero according to e 1 / 2 6{t). The singularity is not an artifact of the simple model, 
and should be observable in real materials. Roundoff may occur due to lifetime effects and sample inhomogeneities. 



G. Ferroelectric resonance 



The dash-dotted line in Fig. ^ represents the absorption spectrum of the mixed- valent system in a bias field. The 
spectrum consists of two parts: a peak at u> = ujq, and a continuum above uj = 2A. We shall show that ujq is the 
ferroelectric resonance frequency. The continuum part of the absorption spectrum is blue-shifted only slightly by the 
bias field. 

Ferroelectric resonance occurs when an electronic ferroelectric, placed in a bias field E z , is acted upon by an 
alternating field £ z of frequency ujq. The ferroelectric resonance frequency ujq is the frequency of the Goldstone mode 
in the bias field E z . The Goldstone mode corresponds to a uniform precession of the pseudo-spins around the z-axis. 
For an ideal system c^o vanishes at zero bias field. For a real system the external bias field E z must be replaced by an 
effective internal field £cff , which may depend on the sample shape, the domain structure, and the crystal anisotropy. 
An example of a contribution to E e g is the depolarization field —N Z P Z of the sample boundary, where N z is the 
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depolarization factor. The effective field remains nonzero in the absence of the external bias field, yielding a finite 
resonance frequency for an unbiased sample. 

Ferroelectric resonance is the electric analogue of ferromagnetic resonance in a ferromagnetic insulator. In the 
magnetic case, the alternating field causes a uniform precession of real spins around the z-direction. The ferromagnetic 
resonance frequency is ujq = jH z , where 7 is the gyromagnetic ratio, and H z is the bias field. In real samples H z 
must be replaced by an effective field H c ff depending on the sample geometry, magnetic domain structure, and crystal 
magnetic anisotropy pj|] . 

Ferroelectric resonance does not occur in displacive ferroelectrics because the order parameter only has a discrete 
symmetry. This means there is no Goldstone mode in the excitation spectrum of a displacive ferroelectric. Ferroelectric 
liquid crystals do have an order parameter with a continuous symmetry. The dielectric response of the Goldstone 
mode has been observed [ 22| in ferroelectric liquid crystals by means of broadband dielectric spectroscopy. 

We now calculate the dependence of the ferroelectric resonance frequency of an ideal system on the external bias 
field E z . For a real system, E z must be replaced by E e g. The ferroelectric resonance frequency is given by the 
equation D(u)q) — 0, where D(u>) is the denominator of Eq. (|28|). For [i z E z small compared to A, an approximate 
solution may be obtained by expanding D(ujq) in a Taylor series around uj = 0: D{uj q ) = D(0) + ±c^L>"(0). The 
linear term vanishes because D{u>) is even in uj. Neglecting terms of order R 2 in D(0) and of order R in D"(0), we 
find 



UJ 



-4m z A I 



A(0) 



4AM 2 (0) + B 2 (0) 



1/2 



(31) 



Substituting the explicit expressions for ^4(0) and B(0) given in the Appendix yields 



UJ 



U 



1/2 



(32) 



where E' c = E c — U is the critical value of E'f . A useful estimate of ujq is the arithmetic mean of the gap 2 A and the 
field energy ix z E z . For fi z = 10 -29 Cm, 2A = 1 meV, and E z between 10 V/cm and 10 4 V/cm, ujq is between 0.01 
mcV and 1 meV, i.e. in the microwave regime. 

The condition D(ujq) = is not sufficient for a peak in the absorption spectrum at uj = uj$. One must also have 
N(ojq) ^ 0, where N(u>) is the numerator of Eq. (p8|). The strength of the pole at ujq is Zq = —■kN{ujq)/D'{u>q). 
Away from half-filling, one finds Z > 0. Exactly at half-filling, B{uj) = 0. As can be seen at once from Eq. (p8|), we 
then have Zq = 0. The zero strength at half-filling is an artifact of our simple model, which has p(e) = p(— e). For a 
realistic d-band there will be a resonance peak at half-filling. 

An important problem in ferromagnetic resonance is to account for the width of the ferromagnetic resonance 
line [ p3| . The line width is due to the relaxation of the uniform precession by spin-spin and spin-lattice interactions. 
Various damping terms may be added to the Bloch equations to describe the relaxation. A possible damping term is 
— (Sk — S^^/t, where r is a phenomenological relaxation time. In the equation of motion for , this leads to the 
replacement of d/dt by d/dt + 1/r. Thus, the simplest approach to damping (the one we have adopted here) is to 
add an imaginary part 8 = 1/r to to. A detailed study of the width of the ferroelectric resonance line is left for future 
research. 



IV. SECOND-HARMONIC GENERATION 



Second-harmonic generation is the generation of an outgoing electromagnetic wave of frequency 2oj from two in- 
coming waves of frequency uj. For incoming waves propagating along ki and k2, the second-harmonic radiation is 
most effectively generated in the phase- matching direction ki +k2. In this direction both energy and momentum are 
conserved. 

The ability of a medium to sustain second-harmonic generation is characterized by the second-harmonic suscep- 
(2) 

tibility tensor x^(2w, uj, uj). Here i,j,k — x,y,z are the Cartesian indices. The polarization Pi(2uj) induced in 

( 2) 

the medium by incoming fields £j{u) and £fc(w) is given by Pi(2uj) = x^-j,(2o;, uj, uj)£j{uj)£k{uj). For a medium with 

(2) (2) 

inversion symmetry, the second-harmonic susceptibility must satisfy X\jk ~ ~X-ljk = ^- Therefore, second-harmonic 
generation cannot occur in media with inversion symmetry. 

In the mixed-valent system, the inversion symmetry is spontaneously broken by the pairing of electronic states of 
opposite parity. This leads to the appearance of a built-in polarization in the ^-direction. If the incoming fields have a 







component along the z-axis, second-harmonic generation can occur. We calculate the second-harmonic susceptibility 

(2) 

Xzzz for incoming fields polarized along the z-axis. For incoming fields polarized at an angle relative to the z-axis, 
the second-harmonic susceptibility is reduced by cos 2 <j>. 



A. Pseudo-spin nutation 

(2) 

We calculate the second-harmonic susceptibility Xzzz from the optical Bloch equations by expanding Sk and 
to second order in £ z . The external "magnetic" field Hk has no components of second or higher order. The equation 

(2) 

of motion for S k is 

SP> - (H k 0) - M k 0) ) x S k 2) + M<?> x S k 0) = (H k 1} - M«) x S«. (33) 
Eq. d33| ) has the same form as Eq. (17), except with a more complicated right hand side. As before, we decompose S k 



(2) (2) (2) (2) 

into its radial, polar, and azimuthal components: S k = S, r k e r + Sg k e# + k e0. To second order, the pseudo-spin 

(2) 

has a nonzero radial component S^ k . This means the motion is no longer a regular precession: the pseudo-spin 
nutates during the precession. (Nutation is the up-and-down motion of the precession axis.) The nutation frequency 

(2) 

is twice the precession frequency. This can be shown as follows: the equation of motion for S^. k is obtained by taking 
the inner product of Eq. (|3^) with e r . Using the vector identity (A x B) ■ C — (C x A) • B and taking the inner 

product 
we find 



product of Eq. (llTl) with e r shows that the equation for S^ 2 k can be written as S^ 2 k = — S^ 1 "* ■ S k ^. Upon integration, 



:( 2 ) „„a c(2) 



2 

The equations of motion for Sg\ and SY'^ are obtained by taking the inner product of Eq. ( p3| ) with eg and e 



( 2) 

respectively. Apart from the more complicated source terms, a new feature occuring in second order is that k is 
now nonzero. This leads to an additional term ^ k , sin(#k — ^k')^rk' on * ne nan( i s fo e °f the equation for S^ k . 
However, since we have already solved for 5^ 2 k in Eq. (|34|), the additional term can be taken over to the right hand 

(2) (2) 

side and treated as an extra source term. Combining all source terms into driving forces F e k and k , the equations 
of motion for Sg 2 ^ and S^ 2 k are 

sfl-2E^l + Mfl = Fll (35) 



with 



^U^sfl-M^ = Ffl (36) 



F e ( % = (2^ z £ z sin k + M)J)S;}, (37) 
Ffl = - (2» z £ z sin k + M$)S$ - ^ £ sin(0 k - . (38) 



k' 



Here = S J^ k , sin(0 k — 0k')^0k'- A very important observation is that since sin 8^ = A/Ek, all source terms 

are proportional to A. The second- harmonic susceptibility vanishes identically when A = 0. 



B. Independent quasiparticles 

We first calculate the second-harmonic susceptibility neglecting the Coulomb interaction between the optically 
excited quasiparticles. This corresponds to setting M^ 2 k = M^ k = on the left hand side of Eqs. ([35]) and (^), and 
Fg k = 2/i z £ z S'^ 1 k sin6' k and F^ — —2fi z £ z Sg 1 ^sm6i !i on the right hand side. The second-harmonic susceptibility is 
Xzzz = P z 2 ^/£ z , where P^ = (JV// z /fl) J^kO^k cos #k + Sy k sin 0k) is the second-order polarization. The second- 
order polarization has a contribution from the radial component Sy k . Solving Eqs. (|3^) and ( |36| ) for Sg k , and using 
Eq. ( |34| ) for 5^ k , we obtain 
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(w + 



(w - 2E k )(2uj - 2E k ) 

2m k ,k« k ,k 2m k,k" k ,k 
2£ k )(2w + 2£; k ) ~ (w- 2£ k )(w + 2£; k ^ 



(39) 



The coherence factor m k , k ' is given by ™ kjk ' = it k v k ' + v k u k r. The physical origin of the coherence factor m k , k < is the 
scattering of a quasiparticle by the second incoming photon. The second photon can either scatter the quasi-electron 
from k to k', or the quasi- hole from — k' to — k. The overall amplitude for the scattering process is 2m kk '. Since 
TO k)k = A/_E k , the second-harmonic susceptibility is directly proportional to A. 



C. Final-state interaction 



For a separable Coulomb potential, an analytic solution for the second-harmonic susceptibility including the final- 



state interaction is possible in principle. However, the large number of driving terms in Eqs. (37) and ( |3q ) presents 
a considerable challenge. We have instead approached the problem numerically. This is done in analogy with the 
classical mechanics treatment of forced oscillations. The azimuthal component S^, tk is the generalized coordinate q k , 
and the polar component Sg tk is the negative of the conjugate momentum p k . Eqs. ( |35|) and ([56]) are the Hamilton 
equations of motion for q k and p k . The matrices T and V are 



U 
N 
U 



T~l = 2E k 5 k . k , - — cos(0 k -fa), (40) 



Vk, k < - 2E k 8 k . k , - — . (41) 

The first step is to find the normal modes of oscillation of the system of pseudo-spins. The normal-mode equations 
are 

Vk,k'^k',„ = (2E n ) 2 T k , k -A kl , n , (42) 

k' k' 

where 2E n is the frequency of normal mode n, and A kn is the amplitude of £<A, k in the normal mode n. Since T 
and V are both real and symmetric, the frequencies are all real and positive. There is one Goldstone mode, whose 
frequency is the ferroelectric resonance frequency u)q. The remaining TV — 1 frequencies form a continuum above the 
energy gap. 

(2) (2) (2) (2) 

The second step is to obtain the forced oscillation of Sg k and k when driven by Fg k and F^ k . This is done by 
solving for the motion in normal coordinates, and then taking linear combinations to obtain the motion in the original 
coordinates. The force driving the normal coordinate £„ has frequency 2oj and amplitude Q n = — ^2 k {A k . n Fg J + 
2iu;A~ 1 k F^l). This causes the normal coordinate to oscillate with frequency 2ui and amplitude £n = Qn/(4£^ — 4w 2 ). 
The original coordinates oscillate with frequency 2oj and amplitudes S\ k = 2j n A kj „£„ and Sg k — 2iiv ^ n A^ n ( n + 
Sk' -^k.k'-Pl k " The complete motion in the original coordinates is (m = 2): 

qim) \- M^jg + rni^A^F^) 

^ k ~ 2^ m 2 w 2 _ AE 2 > ^ 

^ KU^A^,F$- m *.A k , n F$) 

n,k' " 

Here we have used the orthonormality condition S„ tn > — k , A k n T k k i A k ^ n r for the matrix of eigenvectors A k n . 
Eqs. (pl3j ) and (|44|) apply to m-th harmonic generation in general. 

The computational task is summarized as follows: First compute and Sg 1 ^ from Eqs. (43) and (^||) with m = 1 



and Fg 1 ! = 0, F9~! — 2fi z £ z cos9 k . To avoid a singular denominator, u> is given a small positive imaginary part S. 



Then compute S^l from Eq. @, and and F^l from Eqs. @ and @. Finally, compute sfQ from Eq. @) 
with m = 2. 
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The results of the calculation are shown in Fig. |4|. The figure shows the amplitude \x z Jz{2oj,uj,uj)\ of the second- 
harmonic susceptibility as a function of the photon energy oj, for several values of Ef, in zero bias field. The important 
features are: (1) The second-harmonic susceptibility is directly proportional to A. (2) When the /-level lies inside 
the d-band, the second-harmonic conversion efficiency is strongly enhanced at lu = A, and less strongly at u> = 2A. 
The first feature shows that second-harmonic generation can be used to test the validity of the SCMF solution in real 
mixed- valent compounds. The second feature distinguishes the single-quasiparticle treatment of the second-harmonic 
response from the SCMF treatment. The enhancement of the conversion efficiency is due to the final-state interaction 
between the optically excited quasiparticles. 

V. EXISTING EXPERIMENTAL EVIDENCE AND PROPOSED TESTS 

As experimental tests of electronic ferroelectricity in a mixed-valent compound we propose measurements of the 
static dielectric constant, the microwave absorption spectrum, and the second-harmonic susceptibility. As an example, 
consider SmB6. The crystal structure of SmBg has cubic symmetry, with J5q octahedra at the body center, and Sm 
ions at the corners of a conventional bec unit cell with lattice constant a = 4.13A. The crystal has inversion symmetry 
at the bec lattice points. Through measurements of the ionic radius, the valence of the Sm ion is found to be 2.53, 
almost halfway between 2 and 3, so that the /-level lies near the center of the conduction band. 

The measured far-infrared absorption spectrum p|-|Iq] of SmB 6 can be interpreted in accordance with the SCMF 
solution. In Fig. || we compare the mean-field and single-quasiparticle results for the linear susceptibility to experi- 
mental data on SmBg taken from Ref. |u| . The data show an energy gap around 2 A = 4 meV, and a sharp peak at 
threshold. The mean-field theory fits the data very well in the threshold region, whereas the single-quasiparticle theory 
gives a qualitatively wrong behaviour. Away from threshold, discrepancies between mean-field theory and experiment 
occur because of our simple model. Further experimental indication of the validity of the SCMF solution in SmBg is 
provided by the electron tunneling spectrum [p.7| , which can be interpreted by analogy with Giaever tunneling in a 
superconductor. 

Ref. @ also reports a measurement of the static dielectric constant of SmBg at T = AK. The large observed 
value e = 1, 500 provides further support for the theoretical prediction of electronic ferroelectricity. An interesting 
test would be the existence of the ferroelectric resonance in this compound. For the parameter values of Fi g. pL the 
ferroelectric resonance frequency is ll>q — 0.021-^/1^1 meV, where E z is in V/cm. This estimate from Eq. (132) was 
checked by the numerical solution of D(ujq) = and was found to be in agreement to better than 0.5% for fields 
up to 10 3 V/cm. For a reasonably strong applied electric field, the ferroelectric resonance lies easily in the range of 
frequency measured in Ref. p4{ at zero static electric field since the lowest frequency measured there was 1 meV. 

Ref. |24| reports measurements of the dielectric response of mixed-valent Sni3Se4 and Sm.2Se3 from 20 Hz up to 1 
GHz. The huge observed values of the static dielectric constant (e = 30, 000 for Sni3Se4 and e = 4, 000 for Sn^Sea) are 
consistent with electronic ferroelectricity. Given the value 2 A = 140 meV for the energy gap in Sni3Se4, we predict a 
resonance frequency between 0.1 meV (10 GHz) and 10 meV (1 THz). 

Ref. pj| reports transport measurements on mixed-valent TmSeo.45Teo.55 that show evidence for a condensation 
of free carriers into an excitonic insulator ground state. Due to the indirect nature of the energy gap, this material 
might not be suitable for the optical tests of coherence proposed here. 

VI. DISCUSSION OF OTHER GROUND STATES 

We now discuss the possibility of electronic ferroelectricity, ferroelectric resonance, and second-harmonic generation 
for solutions of the FK model other than the SCMF solution. There are two types of solutions we consider here: (1) 
solutions with a classical /-electron distribution, and (2) electronic-polaron solutions. The question we ask is: does 
the theoretical ground state have a built-in coherence between d-electrons and /-holes? This coherence is necessary 
for electronic ferroelectricity and the concomitant optical signatures to occur. The answer is, for the two types of 
solutions: (1) no, and (2) yes. A brief motivation for each answer is given below. 

In Refs. ||,|| the /-occupation fjfi is replaced by a classical variable Wi, where Wi = 1 when site i is occupied 
by an /-electron and Wi = when it is not. With this replacement, the FK model becomes a tight-binding model 
with an on-site potential that can assume two different values U or 0. For a given configuration of /-electrons the 
ground state energy is found by filling the lowest ci-electron levels. The ground configuration is the configuration with 
the lowest ground state energy. The question we ask is: what is the value of A = jj S»=i(^l/< ^ilV')? Since each 
site is either occupied {fjfi — 1) or empty (fjfi — 0), every term in the sum vanishes, and A = no matter what 
configuration of /-electrons we choose. 
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However, we argue that there are degenerate states in the /-occupation representation, linear combinations of which 
are ground states with hnite polarization. It is straightforward to demonstrate this property for the exact solutions 
of a ring of four sites. We construct a ring of four sites with d- and /-states on alternate sites. The Hamiltonian, 
including all the FK terms, has inversion symmetry with respect to the d- and /-sites. There are three eigenstates 
which are even under inversion with respect to a d-site, and three which are odd. The parameters are varied so that 
a degeneracy occurs between two states of odd parity and one of even parity. The ground state in the limit of a 
vanishing electric field is a mixture of even and odd states with a finite polarization. 

In Ref. [|| an electronic-polaron solution of the FK model is presented that may explain the anomalous properties 
of heavy-electron materials such as UPt3, UBei3, and CeCu2Si2- It is proposed that the Coulomb interaction causes 
the /-electrons in these materials to propagate like polarons, with a screening cloud made up of d-electrons. The 
model is shown to account for a large number of thermodynamic and transport properties of heavy-electron materials. 

Here we show that the electronic polaron has a nonzero built-in coherence 

k 

where Gd/(k, t) = — (T T d^(r) f^(0)) is a mixed Green's function with Fourier transform Gd/(k, iu> n ) — 
— VGff(k,iuj n )GddO^,i^n), where Gdd and G// are the d- and /-electron propagators, and uo n are the Matsub- 
ara frequencies. Expressing the Green's functions in the Lehmann representation, and summing over the Matsubara 
frequencies, we find (at T = 0) 

POO 

A = VU ATtcsc{-Ka) I dLuuj- a N + (uj). (46) 

J MO 

Here A and a are the forefactor and the singularity index of the /-electron spectral function (see Ref. ||), N+UjS) is 
the density of d-electron states in the upper branch, and /iq is the Fermi level at T = 0. Using Eq. (8) of Ref. H for 
N + (ui) we obtain 

A=^f^VV, (47) 



where W is the bandwidth and rj the upper-branch threshold. Taking the parameter values a = 0.8, = 1.3 77, and 
U = 20 V from Ref. § we get A = 3.28 W. 



VII. SUMMARY AND CONCLUSION 



In this paper we have investigated the linear and nonlinear optical characteristics of the Falicov-Kimball model. The 
SCMF solution of the periodic model results in the Bose-Einstein condensation of d-f excitons. We found that the 
pairing of d-states of even parity with /-states of odd parity breaks the inversion symmetry of the underlying crystal, 
leading to electronic ferroelectricity. The valence transition is accompanied by a divergence of the static dielectric 
constant at the critical value of the /-level energy. The existence of electronic ferroelectricity in a given mixed-valent 
compound predicates on the dominance of the d-f Coulomb interaction over the hybridization. 

We have calculated both the linear and the second-harmonic susceptibilities of a model mixed-valent system within 
the SCMF approximation. The absorption spectrum of the mixed-valent system, when placed in an additional static 
electric field, consists of a peak at the ferroelectric resonance frequency and a continuum above the energy gap. The 
ferroelectric resonance frequency is proportional to the square root of the effective bias field, which depends, in addition 
to the applied static field, on the sample shape, domain structure, and crystal anisotropy. The continuum absorption 
has a threshold singularity when the /-level lies within the conduction band. The second-harmonic susceptibility 
is directly proportional to the amount of built-in coherence A. The final-state Coulomb interaction enhances the 
second-harmonic conversion efficiency at to = A and u) — 2A. As experimental tests of the electronic ferroelectricity 
in mixed-valent compounds we proposed measurements of the static dielectric constant, the microwave absorption 
spectrum, and the second-harmonic susceptibility. The measured far-infrared absorption spectrum of SmBg was found 
to be consistent with the model calculation. 

We have also discussed the possibility of electronic ferroelectricity, ferroelectric resonance, and second-harmonic 
generation for two other theoretical ground states of the Falicov-Kimball model. The electronic-polaron state does 
have a built-in coherence comparable to the SCMF solution. A ground state with /-occupation as a good quantum 
number has no built-in coherence between d-electrons and /-holes. Such a state on its own is not ferroelectric, does 
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not exhibit ferroelectric resonance, and cannot sustain second-harmonic generation. We have argued that degenerate 
ground states in the /-occupation representation can lead to ferroelectric ground states. The explicit numerical 
demonstration of such ground states is left for the future. 
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APPENDIX 



In this appendix we give the explicit expressions for the functions A(u>) and B(u>) defined in Eqs. ( |26| ) and (|2 
Replacing ± £ k by / dep(e), with p(e) = 6{W - |e|)/(2W), we obtain 



At S U 

A{UJ) = 2W 



W-E 



f 



de 



1 



B{u) 



U 

2W 



■W-Ef 
W~E f 



de- 



Ve 2 + M 2 (io 2 - e 2 
e 



M 2 )' 



-W-Ej We 2 + M 2 (uj 2 -e 2 -M 2 ) 
The integrals can be performed by elementary methods. We find 



U 



4W 2 ujVuj 2 - M 2 



In 



' y/(E f + W) 2 + M 2 Vuj 2 — M 2 + (E f + W)u 
y/{E f + W) 2 + APV^^M^ - {E f + W)uj 



and 



B{w) 



( y/(E f - W) 2 + MVw 2 - M 2 + (E f - W)uj S 
11 \ y/(E f - W) 2 + M 2 y/u 2 - M 2 - [E f - W)u> / 



ln / ^{E f +W) 2 + M 2 +uj \ _ ^ / ^{E f -W) 2 +M 2 +u 

4ir^- 11 ^ y/(E f + W) 2 + M 2 -uj) n \^f{E } -W) 2 +M 2 -u / 



U 



(48) 
(49) 



(50) 



(51) 



The frequency u> is understood to have a positive imaginary part S. The function A(oj) is even in Ef, while B(uj) is 
odd in Ef. This means B(u>) = at half-filling. 

The static susceptibility given by Eq. ( |To| ) involves the functions A(u>) and B(u) evaluated at co = 0. Taking the 
limit to — > in Eqs. (|50"|) and (|5l|) we obtain 



A(0) 



B(0) 



U 



W + E 



f 



W-E 



f 



2W 2 M V y/(W + E f ) 2 + M 2 y/(W-E f ) 2 + M 2 



U 
2W 2 



^/(W - E f ) 2 + M 2 y/(W + E f ) 2 + M 2 



(52) 
(53) 



FIG. 1. /-level occupancy n/ and gap parameter A of the model system as a function of the /-level energy Ef. The 
Coulomb repulsion is U = 3.0W. 



FIG. 2. Static dielectric constant e zz of the model system as a function of the /-level energy Ef. The inset shows the built-in 
polarization as a function of Ef. The Coulomb repulsion is U — 15 meV and the electric-dipole matrix element is fi z = 3.7 10 -29 
Cm. The parameter values were obtained by fitting the absorption spectrum of the model system to experimental data on 
SmB 6 (see Fig. §). 
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FIG. 3. Absorption spectrum ImxlV(w) of the model system as a function of the photon energy u), for various values of the 
/-level energy Ef. The Coulomb repulsion is U — 3.0W. The solid lines show the continuum absorption spectrum in zero bias 
field. The dash-dotted line shows both the ferroelectric resonance peak and the continuum absorption spectrum in a bias field 
of E z = 0.010* /W. 



FIG. 4. Amplitude \)$S x (2w , w,w)| of the second-harmonic susceptibility as a function of the photon energy lj, for various 
ilues of the /-level energy Ef. 
jpulsion is U = 3.0W. The am] 
Fig. I, Nfil/(2QW 2 ) = 82nmV" 



(2) 

values of the /-level energy Ef. The dash-dotted line shows the phase of Xzzz(2u),ui,u>) for Ef = —LOW. The Coulomb 
repulsion is U = 3.0W. The amplitude is given in units of N[i^/(2QW 2 ). For the parameter values given for the solid line in 



FIG. 5. Comparison of the mean-field (solid line) and single-quasiparticle (dash-dotted line) results for the infrared absorp- 
tion spectrum of SmBg to experimental data taken from Ref. [[14J (diamonds). The /-level energy is Ef — 0, the d-bandwidth 
is W = 40 meV, and the Coulomb repulsion is U = 15 meV. The electric-dipole matrix element is fi z = 3.7 10~ 29 Cm for the 
solid line, and fx z = 4.2 10~ 29 Cm for the dash-dotted line. 
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